Fragile to strong crossover coupled to liquid-liquid transition in hydrophobic solutions 
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Using discrete molecular dynamics simulations we study the relation between the thermodynamic 
and diffusive behaviors of a primitive model of aqueous solutions of hydrophobic solutes consisting of 
hard spheres in the Jagla particles solvent, close to the liquid-liquid critical point of the solvent. We 
find that the fragile-to-strong dynamic transition in the diffusive behavior is always coupled to the 
low-density/high-density liquid transition. Above the liquid-liquid critical pressure, the diffusivity 
crossover occurs at the Widom line, the line along which the thermodynamic response functions 
show maxima. Below the liquid-liquid critical pressure, the diffusivity crossover occurs when the 
limit of mechanical stability lines are crossed, as indicated by the hysteresis observed when going 
from high to low temperature and vice versa. These findings show that the strong connection 
between dynamics and thermodynamics found in bulk water persists in hydrophobic solutions for 
concentrations from low to moderate, indicating that experiments measuring the relaxation time in 
aqueous solutions represent a viable route for solving the open questions in the field of supercooled 
water. 

PACS numbers: 64.70. Ja,65.20.-w,66.10.C- 



I. INTRODUCTION 

Over the past decade our understanding of the puz- 
zling behavior of water in the supercooled region has 
greatly expanded [l|. At present, plausible scenarios that 
might fill in the missing pieces of our comprehension have 
been developed through simulations and theoretical mod- 
els. In 1992, using a molecular dynamics simulation of 
the ST2 potential, the existence of a liquid-liquid criti- 
cal point (LLCP) in the supercooled region of water was 
hypothesized [2j. Since that time the LLCP has been 
found in a number of other simulations, see for example 
Refs. • In this conceptual framework the presence of 
the LLCP explains the puzzling anomalies in the behav- 
ior of water. Although the LLCP scenario remains fas- 
cinating, its existence has not been experimentally con- 
firmed, and alternative scenarios have been proposed [6|. 
Experiments are consistent with the existence of a LLCP 
@, Q but proof has been difficult because efforts to reach 
the zone in the phase diagram where the LLCP would 
be located are hampered by nucleation. On the low tem- 
perature side the existence of high density amorphous 
(HDA) and low density amorphous (LDA) states with a 
first-order coexistence line separating them has been ex- 
perimentally proven If a LLCP exists, it is at the 
terminating point of the extension of this line in the su- 
percooled liquid region. [1, [Io| ■ 

In 1996 a strong relationship between the dynamic 
and thermodynamic behaviors of water was reported 



11| . indicating that the well-known singular tempera- 



ture T$ [12J at which extrapolation of various thermo- 
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dynamic and dynamic anomalies diverge, can be identi- 
fied with the crossover temperature Tc of the mode cou- 
pling theory (l3| . This crossover temperature marks the 
fragile-to-strong (FTS) transition that occurs in super- 
cooled liquids approaching the glass transition. Recent 
studies have completed the picture for water by connect- 
ing this FTS transition to the crossing of the "Widom 
line" [1, [ID, [H|. Upon approaching the LLCP, the lines 
of maxima of the different response functions converge on 
this Widom line, which separates water with HDL-like 
features at high temperatures from water with LDL-likc 
features at low temperatures [H, HH . 

Simulations have indicated that the LLCP persists in 
aqueous solutions for concentrations ranging from low to 
moderate both for polar and apolar solutes [l7l - [20| . 
This suggests that focusing on aqueous solutions might 
be a possible experimental strategy in probing the ex- 
istence of a LLCP. In solutions, in fact, the nucleation 
barrier can be tuned and brought below the LLCP 
[l8| . Very recently experiments on water/glycerol mix- 
tures [2l| have shown signs of a liquid-liquid phase tran- 
sition. Nonetheless, thermodynamic measurements are 
more difficult than measurements on dynamics, therefore 
it is timely to prove whether the FTS crossover remains 
coupled to the liquid-liquid transition also in aqueous 
solutions, as this can address the experimental quest on 
the existence of LLCP toward the study of dynamics. We 
note that it is a priori not obvious that the FTS transition 
transfers to aqueous solutions as it is a phenomenology 
present in many glass formers that do not have a LLCP. 
In this paper we use a primitive model able to reproduce 
the phenomenology of water. We present the results of 
discrete molecular dynamics (DMD) simulations on wa- 
ter mixed with hard spheres (HS), in order to determine 
whether the connection between dynamics and thermo- 



2 



dynamics retains its validity in hydrophobic solutions. 

The paper is structured as follows. Simulation details 
are reported in Sec. El The results are presented and 
discussed in Sec. IIIII Concluding remarks are given in 
Sec. EH 



II. SIMULATION DETAILS 

The primitive model used for the results presented in 
this paper is the Jagla model [22| [see Fig. QJa)], an 
isotropic potential with two length scales, a hard-core 
distance a and a soft-core distance &, plus an attractive 
ramp that extends to a cut-off distance c. The parame- 
terization we use, b/a = 1.72, c/a=3 and Ur/Uq — 3.56, 
has been shown to possess bulk water anomalies and a 
liquid- liquid transition [J, |23j . The Jagla model was also 
successful in reproducing the increase of solubility upon 
cooling of non-polar solutes modeled as HS of diameter 



Physical quantities are expressed in reduced units fl9| . 
The total number of particles contained in the cubic sim- 
ulation box is N — 1728. We performed DMD simula- 
tions on the discretized version of the Jagla ramp poten- 
tial, for the bulk and for £hs = 0.10 and xhs = 0.15 
HS(J) mixtures of HS in Jagla particles, at constant N, 
P and T. Pressures and temperatures are controlled 
by Berendsen algorithms. The HS have diameter a and 
the same mass as the Jagla particles. Solute-solute and 
solute-solvent interactions are purely hard-core. The po- 
sition of the LLCP of the bulk Jagla system 0, [H, Hi| and 
of HS(J) mixtures with mole fraction up to xhs — 0.50 
has been previously determined [191 ] in NVT simulations. 
Its position in the P — T plane shifts to higher pressures 
and lower temperatures upon increasing the mole fraction 
of HS. In the Jagla model the slope of the liquid- liquid 
coexistence line is positive, unlike other models for wa- 
ter [l[ . Consequently the HDL is the more ordered phase 
and has an Arrhenius dynamic behavior, while the re- 
verse is true for other models [3, Q11 [23[. Simulations 
with primitive models have the advantage of permitting 
to explore for large systems, a wide range of pressures 
and temperatures, especially in the deep supercooled re- 
gion, where equilibration of long-range, orientational de- 
pendent potentials becomes impossible. In particular, 
the ramp potential with the choice of parameters used 
here is able to reproduce, mutatis mutandis, the complex 
thermodynamic scenario of water commonly derived from 
long-range orientational dependent potentials @, 0, H| . 



III. RESULTS AND DISCUSSION 

In order to verify the connection between dynamic and 
thermodynamic behaviors, we consider the phase dia- 
grams of bulk water and aqueous solutions previously 
determined for the Jagla potential and we study the ther- 
modynamic and diffusive behavior along constant pres- 




FIG. 1. (Color online) (a) Spherically symmetric Jagla ramp 
potential. A snapshot of the xjjs ~ 0.10 system is also shown 
in the top right corner, (b) Schematic liquid-liquid phase di- 
agram for Jagla systems (bulk or solutions) with the LLCP 
and the limit of mechanical stability (LMS) lines. The con- 
stant pressure paths are simulated at P = P c + 0.020 (path 
a); P = P C + 0.015 (/3); P = P C + 0.010 (7); P = P c (6); P = 
P c — 0.010 (e) starting from high temperature; P — P c — 0.010 
(e') starting from low temperature; P — P c — 0.030 (£). 



sure paths [see Fig. [ljb)]. In particular, we consider 
isobaric paths near and above the LLCP, at the LLCP, 
and near and below the LLCP. Paths a, (3, and 7 are 
above the critical pressure and cross the Widom line. 
The path S is at the critical pressure. Paths e and e' 
are below the critical pressure, and cross the HDL and 
LDL limit of mechanical stability (LMS) lines depending 
on the direction. In particular, along path e the system 
goes from LDL (high T) to HDL (low T) when crossing 
the LDL LMS while along path e' it goes from HDL to 
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FIG. 2. (Color online) Isobaric specific heat Cp as a function of temperature (a,b) of bulk Jagla particles; (c,d) of the xhs = 0.15 
HS(J) solution. Top panels are for paths with P > P c , bottom panels for paths with P < P c . Lines are cubic splines to the 
points. 



LDL upon crossing the HDL LMS. Path £, also below 
the critical pressure, does not cross the LDL LMS line 
and therefore the system remains in the LDL all along 
the path. All paths, except path e', have been performed 
starting from a configuration at high temperature and 
progressively equilibrating the system at lower tempera- 
tures. For path e' the starting configuration was taken at 
low temperature and the system was progressively equi- 
librated at increasing temperatures in order to cross also 
the HDL LMS. 

Figure [2] shows the behavior of the isobaric specific 
heat C P = (1/N) ■ (dH/dT) P for bulk water and for the 
%hs = 0.15 HS(J) solution calculated for all the paths 
shown in Fig. [ljb). The results for the xns = 0.10 
HS(J) solution follow a similar trend and are not shown. 
Both in bulk water and in the solutions, the points above 
the critical pressure at which the specific heat displays 
a maximum define the Widom line emanating from the 
LLCP. Consistent with the positive slope of the liquid- 
liquid coexistence line for the Jagla potential 

mm, the 

specific heat maximum moves to a higher temperature 
as the pressure increases. Its height decreases moving 
away from the LLCP. For the paths along the critical 
pressure, the temperature of the maximum falls close to 
the one estimated in NVT simulations [TH, [23| ■ We note 
that paths e and e' clearly show the hysteresis expected 



for a first-order transition in a finite size system with 
the Cp maxima occurring at different temperatures, one 
corresponding to the LDL LMS for path e and one cor- 
responding to the HDL LMS for path d . Along path £, 
which does not reach the LDL LMS, no peak is observed. 

We now show the diffusive behavior of Jagla parti- 
cles in solution along the same constant pressure paths 
and compare it with the bulk Jagla results. Figure [3] 
shows the diffusion coefficient as a function of inverse 
temperature for bulk Jagla particles and xhs = 0.10 and 
0.15 HS(J) solutions. For all three systems, above the 
LLCP the diffusion coefficient shows a crossover as it 
crosses the Widom line as determined from the Cp max- 
ima. Along path £, that does not cross the LDL LMS 
line, the temperature dependence of D is clearly non- 
Arrhenius, indicating that the LDL liquid can be classi- 
fied as fragile. Taking advantage of this particular poten- 
tial we could push this path to very low temperatures and 
thus observe a clear fragile behavior. In Fig. [3] we also 
show the fit to the Vogel-Fulcher-Tamman, VFT, func- 
tion D = Dq exp[— Bj (T — To)] characterizing the fragile 
behavior of a glass former. 

For all the paths crossing the Widom line, on the high 
temperature side of this line we find a LDL-like liquid, 
and the diffusive behavior is that of a fragile liquid, when 
compared to path £ which is always in the LDL fragile 
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FIG. 3. (Color online) Diffusion coefficient D as a function of inverse temperature (a,b) of bulk Jagia particles; (c,d) of the 
^HS = 0-10 HS(J) solution; (e,f) of the xjjs = 0.15 HS(J) solution. Top panels are for paths with P > P c , bottom panels are 
for paths with P < P c . The solid curve is a VFT fit D = Do exp[— B/(T — To)] for paths The dashed straight line is an 
example of Arrhenius fit D — Do exp(— Ea/ ksT) for paths e. The parameters of the fits are Do — 0.4, B ~ 0.25, To ~ 0.18 for 
VFT and Do ~ 10 s , Ea — 6 for Arrhenius, for all systems. 



environment. On the low temperature side of the Widom 
line the diffusive behavior shows a steep straight slope 
typical of a strong liquid with D following the Arrhenius 
law, D = D cxp(-E A /k B T). 

Along paths S, at the critical pressure the diffusiv- 
ity jumps from high-temperature fragile behavior to low- 
temperature strong behavior close to the estimated tem- 
perature of the LLCP. An analogous jump is found when 
crossing the LDL LMS line along paths e. We note that 
the jump appears only when passing the LLCP or a LMS 
line. In analogy with the hysteresis behavior observed 
for the isobaric specific heat (see Fig. a jump is also 
observed when heating up the systems along paths e', 
but now it occurs at a temperature corresponding to the 
HDL LMS line. This indicates that the hysteresis caused 
by the first-order transition observed when studying the 
thermodynamic behavior in a finite system is observed 
for the diffusive behavior as well. This jump in the dif- 
fusion coefficient is a finding that can be useful for the 
experimental determination of the liquid-liquid coexis- 
tence. These findings confirm the coupling between the 
dynamic and the thermodynamic behavior observed in 
bulk water (P > P c corresponds to P < P c and vice 
versa in other models [27[ for the reversed slope of the 
coexistence line) and extends its validity to hydrophobic 
solutions. 

Figure H shows the LLCP, HDL, and LDL LMS lines, 
the points at which the isobaric specific heat Cp and 
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FIG. 4. (Color online) Liquid-liquid phase diagram for bulk 
Jagla particles, xhs = 0.10 and xus = 0.15 HS(J) solutions. 
Along with the LLCP (bold large circles), LDL LMS (triangles 
up) and HDL LMS (triangles down), are shown the maxima 
of the isobaric specific heat Cp (small circles), the maxima of 
the coefficient of thermal expansion ap (+) and the crossover 
points of the diffusion coefficient D (x). 
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coefficient of thermai expansion ap maxima occur and 
the points at which a crossover in the diffusive behav- 
ior is found. The coefficient of thermal expansion ap — 
(1/V) • (dV/dT)p has been calculated for all our systems 
and for all paths. It displays a behavior qualitatively 
equivalent to that of the isobaric specific heat. Above the 
LLCP along the Widom line Cp and ap exhibit maxima, 
and the temperature at which the diffusivity crossover oc- 
curs, estimated from the derivative of the logarithm of D 
with respect to 1/T, coincides with Cp and ap tempera- 
ture maxima and thus with the Widom line. The maxima 
in thermodynamic response functions and the diffusivity 
crossover at the critical pressure are also found close to 
the LLCP temperature, estimated in NVT simulations 
[HI, [23[ ■ Below the LLCP and along paths e and e' the 
points where Cp and ap and the diffusivity crossover oc- 
curs also coincide, and they correspond to the crossing of 
the LDL LMS line and the HDL LMS line, respectively. 
This shows that the hysteresis expected for the thermo- 
dynamic quantities calculated along paths e and e' due 
to the first-order coexistence line between LDL and HDL 
and observed only at the LMS lines due to the metasta- 
bility of the phases, is also maintained in the diffusive 
behavior. 



IV. CONCLUSIONS 

Our results suggest that the FTS dynamic transition 
reported for water in different environments in the lit- 
erature [1 [H O GJ, H HI] occurs at the same time 
as the thermodynamic liquid-liquid transition. They are, 



in other words, two sides of the same coin and, when 
the isobaric paths do not cross the LMS or the Widom 
line, there is no crossover and the dynamic behavior re- 
mains strong or fragile. We note that this connection be- 
tween dynamic and thermodynamic behaviors might be 
a general feature of network- forming liquids, e.g., a simi- 
lar picture is found in silica [HI . Why the FTS dynamic 
crossover is also commonly found in glass formers [30l | 
when it is not associated to a liquid-liquid transition re- 
mains an open question. Our results show that for mod- 
erately concentrated solutions the changes induced by the 
solutes are continuous and the FTS crossover remains 
coupled to the liquid-liquid transition paving the way 
to the experimental exploration of aqueous solutions as 
tools to understand the mysteries of water. In particular 
light scattering and neutron scattering experiments (3lj 
can be used to detect the liquid-liquid transition in aque- 
ous mixtures. Particularly important would be to see the 
jump in relaxation times below the LLCP as it is directly 
connected to the liquid-liquid transition. 
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